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We present a three-band Hubbard Hamiltonian and the associated Cu K-edge RIXS spectra 
for electron and hole doped cuprates over a wide range of energy and momentum transfers. By 
comparing computed spectra for the unfilled case with the corresponding results for 15% electron 
or hole doping at two different values of the effective Hubbard parameter U, generic signatures of 
the collapse of the magnetic gap and the characteristic momentum dependencies and evolution of 
the spectra with doping are identified. Available RIXS data support the gap collapse scenario for 
electron doped cuprates, but the situation in hole doped systems is found to be less clear. 

PACS numbers: 



In a resonant inelastic x-ray scattering (RIXS) exper- 
iment, the incident x-ray photon is scattered after in- 
ducing electronic transitions in the vicinity of the Fermi 
energy (Ef), as in Raman scattering)! Since no charged 
particles enter or leave the sample, the technique is at- 
tractive as a genuine bulk probe and a new spectrocopic 
tool for getting a handle on the nature of both filled and 
unfilled electronic states within a few eV's of Ep- Al- 
though filled states are also accessible via angle-resolved 
photoemission spectroscopy (ARPES), the unique abil- 
ity of RIXS to access unfilled states makes it ideal for 
probing how the magnetic gap in insulating half-filled 
cuprates evolves with doping. The momentum-space evo- 
lution of this process and the possible collapse of the 
magnetic gap at optimal doping are fundamental pieces 
of information necessary for identifying different routes 
from the insulating to the superconducting state and the 
mechanism which underlies the remarkable phenomenon 
of high-temperature superconductivity. 

This article discusses a three band Hubbard Hamilto- 
nian and the associated Cu K-edge RIXS spectra for en- 
ergy transfer ui extending to 8 eV and momentum trans- 
fer q along the (tt, it) — > T — > (ir, 0) symmetry line. Our 
purpose is to obtain a systematic scheme for interpreting 
spectral features and their evolution with electron and 
hole doping in the cuprates. The specific parameters are 
chosen to describe Nd2- a: Ce a ;Cu04 (NCCO), but the re- 
sults are expected to be representative of the cuprates 
more generally. The electronic structure is treated at the 
Hartree-Fock level with variable U to simulate the col- 
lapse of the magnetic gap. The solution at half- filling for 
U = 7.2 cV is found to be antiferromagnetically (AFM) 
ordered and to display a 'charge transfer' gap of « 2 eV 
between the Cu d and O-p orbitals. At fixed [7=7.2 eV, 
the effect of electron doping is to move most RIXS tran- 
sitions to lower energies (in relation to half-filling) , while 
hole doping introduces new 'bands' of transitions in the 
spectrum, although both electron and hole doping yield 
additional weak transitions at low energy transfers in the 
0-2 eV range. For [7=5.7 eV, on the other hand, where 
the AFM gap collapses near optimal electron or hole dop- 
ing, strong features in the RIXS spectra are predicted 



with characteristic q-dependencies at low energies. In 
confronting our predictions with available experimental 
resulta 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 : 9 , we find that the RIXS data on NCCO 
support the scenario of a collapsed magnetic gap at op- 
timal doping, while the situation with the hole doped 
systems is less clear. Our study suggests regions of uj — q 
space where further experimental work would be valu- 
able. 

As general background, we note that Cu _ftT-edge RIXS 
is an indirect probe of the magnetic gap, since the x-ray 
excites a Is — > Ap core electron, leading to a shake-up 
of the near-Fermi level states associated with screening 
the core hole. The relative merits of indirect and direct 
RIXS in the cuprates have been debatecUSaii* 1 ^, where 
the latter refers to L-edge RIXS and the excitation of an 
electron from the Cu 2p directly into the Cu 3d states. 
The polarization selection rules in RIXS are equivalent to 
those in Raman. Screening excitations in K-edge RIXS 
satisfy A^=0,2 selection rule, so that the RIXS gap can 
be different from the optical gap (AZ=l)Ai. 

Most treatments of RIXS spectra are based on clus- 
ter model£ 3 - 4-5 j2ii&i!ii£, which are intrinsically limited in 
describing solid-state band dispersions, although some 
studies have involved single band, infinite dimensional 
models 16 . We consider a 3-band Hubbard model Hamil- 
tonian with Cu d x 2_ y 2 and two O p — a orbitals: 

H = ) (A d]dj + Un^riji) + ^ [7 p n lT n a 
3 i 

+ tcuoid^ + (c.c.)} + too\pipi> + (cc.)].(l) 

<i,j> 

Here Ao is the bare on-site energy of the Cu orbital with 
respect to that of the O orbitals (taken to define the en- 
ergy zero), tcuO and too are the Cu-0 and 0-0 nearest 
neighbor (NN) hopping parameters and U [U p ] is the Cu 
[O] on-site Coulomb repulsion. The single [double] angu- 
lar brackets on the summations denote that these sums 
are restricted to the Cu-0 [O-O] NN terms. The spe- 
cific values of the parameters have been obtained from an 
analysisi2ii2*iS of the doping-dependent ARPES spectra 
of NCCO: t C uO = 0.8 eV, t 00 = -0.4 eV, A = -0.755 
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FIG. 1: (color) (a) Density of states (DOS) and its Cu (red) 
and O (blue) components at half-filling for U = 7.2 eV. (b-f) 
Maps of RIXS intensity spectra in the oj-q space approxi- 
mated by WK-JDOS for: (b) half-filling; (c) 15% holes, U = 
7.2 eV; (d) 15% electrons, U = 7.2 eV; (e) 15% holes, U = 
5.7 eV; and (f) 15% electrons, U = 5.7 eV. Dotted lines in (d) 
and (f) highlight prominent low energy features. 



eV, Up = 5.0 eV, and U = 7.2 eV for 1 = 0. £7 is 
found to decrease with doping with a value of 5.7 eV at 
optimal doping, while the other parameters are doping 
independent. Accordingly, computations with U=5. 7eV 
as well as U =7.2 eV are presented. Since similar parame- 
ter values describe LSCOiS, we expect our computations 
to capture generic features of the RIXS spectra of both 
electron and hole doped cuprates. 

As a prelude to RIXS computations, we solve the elec- 
tronic structure at the level of the Hartree-Fock approx- 
imation. Note that the treatment of terms involving the 
Hubbard parameters U and U p , renormalize Ao to an ef- 
fective parameter A = Ao + Una/2 — U p n p /2, where rid 
and md = (n<ff — n^)/2 are the average electron density 
and magnetization on Cu, with similar definitions for the 
oxygens. At half-filling, the self-consistent magnetization 
is found to be nid = 0.295 on Cu-sites, while the O-sites 
are unpolarized. The computed doping dependence of 
nid is similar in 3-bancU^ and one-band^ models. 

For orientation, Fig. ^a) shows the density of states 
(DOS) for £7=7.2 eV at half-filling (a; = 0). The AFM 
ordering yields a total of six bands marked A-F. At half 
filling Ep lies in the gap between the two highest bands 
— the 'upper Hubbard band' A at ~ 3 eV and the 'charge 
transfer' or 'Zhang-Rice' band B at ~ 1 eV. Note that 



the splitting between bands A and F, which are dom- 
inated by the Cu contribution (red), is the main Mott 
gap controlled by U, while the A-B splitting is related to 
the product Urrid, which is the AFM gap 20 . 

The K-shell RIXS cross-section for the Cu Is — > Ap 
core level excitation can be written as2i*2£ 

W(q,w,w0 = (2tt) 3 A|u;(w,^)| 2 
x S(w + Ej (k) - E r (k + q}) nj (k) [1 - n r (k + q)] 
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Here D(u), k) = id + ei s — e 4p (k) + iT ls . The Cu Is band 
is assumed dispersionless at energy ei s , the Cu Ap band 
£4p(ki) is modelled by a 2D tight binding band with NN 
hopping, and the near-FS bands are described by ener- 
gies Ej(k), occupation numbers Uj(k), and eigenvectors 
Xjp of spin er and orbital £ which include a magnetic 
structure factor, oji [u>f] and qj [qf] are the energy and 
wave number of the incident [scattered] photon, with en- 
ergy and momentum transfer given by, to = uji — uif and 
q = qi qf . Polarization effects arise via the Cu Is — Ap 
transition factor 7, which for simplicity is taken as a con- 
stant. The core hole broadening is treated phenomeno- 
logically by the damping parameter r ls with a value of 
0.8 eV^*2i. N is the total number of Cu atoms. Vd is the 
effective potential for creating an on-site Cu 3d electron- 
hole pair. Although the predominant excitations are ex- 
pected to be on-site Cu, we include the possibility of NN 
O excitations via a± = Vi/Vd, where V\ is the potential 
for O 2p excitations, as well as second NN Cu excitations 
via a2 = Vz/Vd- We assume small values of ai = 0.1, 
U2 = 0.05 in this study. is a vector from the core hole 
to the atom on which the electron-hole pair is excited. 
Finally, K a ^ a / is a vertex correction describing excitonic 
coupling of the electron-hole pair; we set h. a ^' — 
which is generally a good approximation^!. 

If the matrix element (ME) w(u>,0Ji) of Eq. [21 is ap- 
proximated as a constant, then Eq. [21 becomes a rela- 
tively simple weighted k-resolved joint density of states 
(WK-JDOS) involving bands j and f where the weights 
involve only Cu— >Cu or O— >0 transitions. More gener- 
ally, at a fixed photon energy ojj, the factor \w(ui, tOi)\ 2 
provides an w-dependent modulation of the WK-JDOS. 
Computations with and without the RIXS cross-section 
approximated as WK-JDOS will be considered. 

Fig. 1 shows the RIXS cross-section approximated as 
WK-JDOS for a number of different cases to allow the 
unfolding of effects of electron or hole doping and those 
due to the collapse of the AFM gap. We focus first on 
the RIXS 'bands' 23 in (b) for the undoped system and 
their relationship to the six DOS 'bands' A-F in (a). Here 
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each of the five filled bands B-F can be excited into the 
unfilled band A, yielding five RIXS bands in (b). The 
lowest RIXS band comes from B — ► A transitions and lies 
around 2 eV, while the highest band around 6 eV in (b) 
comes from F — > A transitions. The overlapping complex 
of bands around 3.5 eV is related to transitions from the 
O-related DOS bands C-E. With adequate resolution and 
statistics, RIXS should be able to map the RIXS bands 
of the sort shown in (b) and thus adduce the nature of 
the ne&r-Ep electronic structure in considerable detail. 

We next consider the effect of electron doping by com- 
paring (b) and (d), where (d) refers to 15% electron dop- 
ing at [7=7.2 eV. Electron doping adds electrons to the 
bottom of the DOS band A, causing the portion of band 
A available for transitions to become narrower, so that 
transitions from bands B-E — > A will display a smaller 
spread in energy. At a fixed value of wid this will induce a 
shift in the RIXS bands to higher energies, but this effect 
is counteracted by a reduction of md with doping (even 
when [7 is held fixed). For example, the B — > A RIXS 
band near T in (d) finally lies at a lower energy of rj 0.2 
eV than in (b). Also, the added electrons at the bottom 
of A can be excited intraband within A and yield a new 
RIXS band at low energies below 1 eV in (d), especially 
along the T - M(tt, 0) line. 

The effect of hole doping is fundamentally different in 
that holes remove states from the top of band B. Since 
bands A and C-F remain unchanged, transitions from C- 
F — > A will see little change, but the narrowing of band 
B due to added holes would move B — > A transitions to 
higher energies (if m<j remained constant)^* Moreover, 
the added holes allow the formation of four new RIXS 
bands associated with excitations from C-F to the top of 
B as seen, for example, by comparing (b) and (c) at the 
middle of the T — M line around 1.4 eV and at 4.3 eV. 
New intraband transitions within B also become possi- 
ble, and much like the case of electron doping discussed 
above, these appear at low energies below 1 eV. How- 
ever, the k-dependencies of such low energy RIXS bands 
for the hole and electron doping are seen by comparing 
(c) and (d) to differ substantially, reflecting the differing 
k-space locations of the respective band edges. 

While md decreases with doping at fixed [7, the result- 
ing magnetic gap is still too large to explain the ARPES 
spectra of NCCO, suggesting that the mean field theory 
does not properly account for screening of U. Reducing 
[7 to 5.7 eV leads to a near collapse of the band splitting 
with m d = 0.063, in good agreement with NCCO ARPES 
data. [For hole doping, [7=5.7 eV yields m rf =0.12.] The 
effects on the RIXS spectrum can be seen by comparing 
Fig. 1(d) with (f) for electron doping and Fig. 1(c) with 
(e) for hole doping. The most prominent signature of the 
collapse of the AFM gap is that the lowest RIXS band is 
shifted down to nearly zero energy at F. 

Fig. illustrates the nature of the factor w(uj,LUi) of 
Eq. 3 and its effect on the spectra, (a) shows the behavior 
of \w\ 2 as a function of <jj for several values of u. [Note 
that Eq. 3 for w does not depend on q.] The peaks in 




FIG. 2: (color) (a) RIXS intensity factor \w\ 2 as a function 
of LOi for several values of u. u>io is the Is-Ap threshold energy. 
Green shaded area schematically shows the DOS for the Cu 
4p band, (b) RIXS intensity map at LUi = uiio for half-filling 
at [7=7.2 eV. 
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FIG. 3: (color) (a) Experimental NCCO RIXS data from Ref. 
.St, Prominent low energy features marked by dotted lines in 
Figs. 1(d) and (f) for 15% electron doping are reproduced for 
[7 = 7.2 eV (+'s) or 5.7 eV (x's). (b) Spectrum of Fig. 1(f) 
on a logarithmic scale to highlight low intensity features, (c) 
Experimental RIXS features in Hgl201 (red dots) from Ref. :9 
superposed on the RIXS map of Fig. 1(c) for 15% holes at U 
= 7.2 eV. 



\w | 2 are seen to be tied to the Van Hove singularities in 
the underlying Ap DOS (green shaded area). In (b), the 
full RIXS spectrum of Eqs. 1 and 2 for x = 0, including 
the effect of w, is shown for a particular photon energy. 
A comparison of Figs. 2(b) and 1(b) indicates that the 
modulation of the spectra due to \w\ 2 is relatively weak. 

Interestingly, the AFM order shows up in the U =7.2eV 
spectra in Figs. l(b)-(d) and Fig. 2(b) as an approximate 
symmetry in the RIXS bands about the (ir/2, tt/2) point 
along the M~r line. However, for [7=5.7eV, the reduced 
magnetism leads to a loss of this symmetry as seen from 
Figs. 1(e) and (f). Note that the underlying bands in all 
cases possess the symmetry of the AFM order—, but that 
this symmetry is in general not obeyed by various spec- 
troscopic measurements [e.g., the 'shadow bands' seen in 
ARPES]. 

We comment now on the implications of some of the 
available experimental RIXS data in the light of our com- 
putations. In this connection, the RIXS data of Ref. @ 
from optimally doped NCCO is reproduced in Fig. 3(a), 
where the lowest RIXS bands highlighted in Figs. 1(d) 
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and 1(f) are superposed as +'s and x's, respectively. In 
comparing the theoretical spectra of Figs. 1(d) and 1(f) 
in the low energy regime of up to w 3 eV with the re- 
sults of Fig. 3(a), it is clear that the experiments are in 
substantially better accord with the results of Fig. 1(f) 
than 1(d). In particular, the data of Fig. 3(a) display 
a larger dispersion and presence of intensity extending 
to quite low energies along (ir, it) — > T, while along the 
r — > (%, 0) line the data show a smaller dispersion and 
are concentrated in a narrower window of energies, rather 
like the results of Fig. 1(f). In sharp contrast, Fig. 1(d) 
shows little intensity in the 0-1.5 eV range. Since the 
computations of Fig. 1(f) for U=5.7 eV involve a col- 
lapsed AFM gap, while those of Fig. 1(d) for U=7.2 eV 
do not, we are lead to adduce that the RIXS data on 
NCCO support an AFM gap collapse scenario. 

Fig. 3(a) also shows a feature around 2 eV at T. 
The intensity of this feature, assigned to interband tran- 
sitions, is enhanced in the representation of Fig. 3(a) 
where data are normalized to maximum at each fc-point 
after subtracting high energy contributions. In reason- 
able accord, a weak feature around 2.5 eV at T is seen 
in our computations in Fig. 3(b), which is due to inter- 
band transitions from the top of the C-E band complex 
to the upper magnetic band^2&. However, we should 
keep in mind that a quantitative description of the RIXS 
spectra with increasing energy transfer will require the 
inclusion of additional Cu-0 and O bands beyond the 
present three-band model. For this reason, other high 
energy features present in the spectra of Ref. |£| and in 
the theoretical spectra of Fig. 3(b) are not analyzed fur- 
ther here. 

The situation in hole-doped cuprates is less clear. Ear- 
lier RIXS experiments^ reported a few, strongly dis- 
persing peaks, but recently Lu, et al.fS find several addi- 
tional weakly dispersing bands both in optimally doped 
HgBa2Cu04+,5 (Hgl201) and in a reanalysis of the LSCO 
data. Fig. 3(b) shows the measured peak positions in 
Hgl201 overlayed on the computed RIXS spectrum of 
Fig. 1(c) for [7=7. 2eV. The placement in energy of the 



observed peaks (red dots) is in qualitative accord, but it 
is difficult to make any definitive assigments of spectral 
features in view of the rather sparse set of q- values sam- 
pled in the experiment. With reference to Fig. 1(e) for 
small U , the most discriminating test of AFM gap col- 
lapse is the position of the prominent dispersing feature 
at low energies along the T — ► (tt/2,tt/2) line. While a 
recent STM studySi hints that the collapse of the AFM 
gap is far from complete in at least some compounds, op- 
tical spectra of hole-doped cuprates find large changes in 
the lowest gap with doping 2 ^. Some differences between 
the optical and RIXS spectra may arise due to presence 
of the core hole in the latter case. Nanoscale phase sep- 
aration may also play a role. The appearance of new 
spectral features with hole doping (e.g. at 4.3 eV and 
at 1.4 eV at (tt/2,tt/2) in Fig. 1(d)), delineated above 
in connection with the discussion of Fig. 1, would be 
tell-tale signatures of the uniformly doped state. 

In conclusion, our study provides a systematic scheme 
for interpreting features in the RIXS spectra of the elec- 
tron as well as hole doped cuprates, although the issue 
of nanoscale phase separation needs to be addressed in 
the latter case. In this way, signatures of collapse of the 
AFM gap in the RIXS spectra can be identified and their 
evolution with doping and momentum transfer vector at 
low energies can be connected to the behavior of the spec- 
tra at high energies to develop a robust understanding of 
the filled and unfilled states within a few eV's around the 
Fermi energy. 
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